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We report on a search for the lepton-flavor-violating processes e + e~ —> fi + r~ and e + e~ — > 
e + r~ . The data sample corresponds to an integrated luminosity of 211 fb -1 recorded by the 
BABAR experiment at the SLAC PEP-II asymmetric-energy B Factory at a center-of-mass energy of 
^fs — 10.58 GeV. We find no evidence for a signal and set the 90% confidence level upper limits on 
the cross sections to be <r MT < 3.8 fb and a eT < 9.2 fb. The ratio of the cross sections with respect 
to the dimuon cross section are measured to be <7 MT /o> M < 3.4 x 10~ 6 and a eT /(7 MM < 8.2 x 10~ 6 . 

PACS numbers: 13.35.Dx, 14.60.Fg, 11.30.Hv 



Within the Standard Model (SM), the fermion mass 
matrices and the mechanism of electroweak symmetry 
breaking remain unexplained. Lepton-flavor is not a con- 
served quantity protected by an established gauge princi- 
ple. Extensions to the SM which include our knowledge 
of neutrino masses and mixing [l[ predict lepton-flavor- 
violation (LFV) at a level many orders of magnitude be- 
low the current experimental sensitivity 0]. 

Searches for LFV have primarily concentrated on the 
decay of the lepton. Limits in a number of muon decay 
channels have reached the 10~ n — 10 -12 level [3| while 
recent measurements of LFV in tau decays have placed 
limits on the branching fractions B(t^ — > /i ± 7) < 6.8 x 
lfT 8 and B(r ± -> e ± 1 ) < 1.1 x 1(T 7 [4] at the 90% 
confidence level (CL). 

There are theories that suggest lepton-flavor can be 
conserved in lepton decay but still be present in produc- 
tion. Some of these models allow for channels such as 
e + e~ — > /i + t~ and e + e~ — > e + r~ through the Q 2 evolu- 
tion of the off-diagonal elements of the fermion mass ma- 
trices §| . Experimental limits on LFV in production are 
considerably weaker than for decay. At center-of-mass 
(CM) energies, y/s — 29 GeV, there are limits on the 
cross section ratios ow/cr^ < 6.1 x 10~ 3 and ovr/cr^ < 



1.8 x 1(T 3 (95% CL) [6[; at V* = 92 GeV, where Z° ex- 
change dominates, B(Z° -> fiT,er) < 0(1) x 1(T 5 (95% 
CL) [7(. The best limits from searches at LEP energies 
above the Z° peak are <t mt < 64 fb and a eT < 78 fb (95% 
CL) No equivalent measurements exist at the lower 
energies accessible by the BABAR detector. 

We present results on two modes of the process 
e + e~ — > 1 + t~ , where l + is an electron or muon and 



the t decays either to n n 



or 7r v T , using data 



recorded by the BABAR detector at the SLAC PEP- 
II asymmetric-energy e + e~ storage rings. Inclusion of 
the charge-conjugate reaction e + e~ — * 1~t + is assumed 
throughout this paper. The data sample corresponds to 
an integrated luminosity of L= 211 fb -1 recorded at a 
CM energy of ^/s = 10.58 GeV. 

The BABAR detector is described in detail in Ref. 
Charged particles are reconstructed as tracks with a 5- 
layer silicon vertex tracker and a 40-layer drift chamber 
(DCH) inside a 1.5 T solenoidal magnet. An electromag- 
netic calorimeter (EMC) is used to identify electrons and 
photons. A ring-imaging Cherenkov detector (DIRC) is 
used to identify charged hadrons and provides additional 
electron identification information. Muons are identified 



by an instrumented magnetic- flux return (IFR) . 

Monte Carlo (MC) simulation is used to evaluate the 
background contamination and selection efficiency. The 
simulated backgrounds are also used to cross-check the 
selection optimization procedure and for studies of sys- 
tematic effects; however, the final background yield esti- 
mation relies solely on data. The signal e + e~ — > Z~t + 
channels are simulated using EvtGen [loj ] in which pho- 
ton radiation is handled by the PHOTOS package 11 1 
to an accuracy better than 1%. The background r-pair 
events are simulated using the KK2F MC generator 1121 1. 
The r decays are modeled with Tauola [l3( according to 
measured rates with the decay t~ — > 7r~7r + 7r _ ^ T assum- 
ing an intermediate aj~(1260) axial- vector state (3. IT3|. 
We also generate light quark continuum events (e + e~ — > 
qq, q = u,d,s), charm, dimuon, Bhabhas, BB and two- 
photon events The detector response is sim- 
ulated with GEANT4 [l6| and all simulated events are 
reconstructed in the same manner as data. 

The signature of the signal process in the CM frame 
is an isolated high-momentum muon or electron recoiling 
against either one or three charged pions and no neutral 
particles. The reconstructed mass of the missing neu- 
trino should be consistent with a massless particle and 
the invariant mass of the recoiling pions and neutrino 
consistent with that of the r. 

We search for events with zero total charge and ei- 
ther two or four well-measured charged tracks originat- 
ing from the e + e~ interaction region. All charged tracks 
must be isolated from neutral energy deposits in the 
EMC and be within the acceptance of the EMC, DIRC 
and IFR to ensure good particle identification. One track 
must be identified as either an electron or muon with a 
CM momentum greater than 4.68GeV/c and no other 
track identified as a kaon or lepton. The electron mo- 
mentum is corrected for energy loss from Bremsstrahlung 
emission by including in the electron momentum the en- 
ergies of isolated calorimeter deposits consistent with a 
photon within a cone of radius 0.1 rad around the initial 
track momentum vector. 

In the CM system, the event topology must be con- 
sistent with an e + //i + recoiling against the remaining 



tracks. We calculate the thrust axis [17| using all the 
charged and neutral deposits in the event and define 
two hemispheres with respect to the plane normal to the 
thrust axis and require that the e + //i + and the other 
tracks to be in separate hemispheres. 
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The t has a fixed CM energy and momentum: 



El = 



r s (M T 2 -Mf) 
2 + : 



P * |= y/E?-M? (1) 



where M r and M; are the masses of the r and e + //i + , 
respectively 0. We define the direction of the r as 
opposite to that of the e + //i + and assign it the mo- 
mentum from Equation [TJ The CM four-momentum of 
the missing neutrino from the r decay, p*, is defined as 
p* — p* , where p* is the sum of the CM four- momenta 
of the pions. The reconstructed r mass is defined to be 
m T = V / (E*+ | p* |) 2 - | p* | 2 where E* is the CM en- 
ergy of the pions. 

Events are rejected if the quantity AE, the difference 
between the e + //i + CM energy and y/s/2, is less than 
—0.5 GeV or greater than 0.2 GeV. True signal events 
will have AE ~ —0.15 GeV while e + e~ — * or 
e + e~ — ► e + e~ events will peak at zero and e + e~ — * t + t~ 
background events have large negative AE. The AE res- 
olution is approximately 50 MeV. Events with converted 
photons are also rejected, where a converted photon is de- 
fined to be a pair of oppositely charged tracks assumed 
to have the electron mass and coming from a vertex with 
a combined mass less than 150 MeV/ c 2 . 

We use a number of kinematic variables to suppress 
backgrounds. The missing event energy in the CM frame, 
E^ iss7 defined as the difference between y/s and the sum 
of the charged track energies, is distributed uniformly for 
signal but peaks at zero or near \fs/2 for the most im- 
portant backgrounds. The missing mass squared, w^ iss , 
should be consistent with zero. A requirement on the 



maximum neutral energy cluster in the detector, 
eliminates events with neutral pions or photons 



E. 



18] 



7> 

A 

requirement on the angle in the CM between the direc- 
tion of the neutrino and the beam axis in the e~ beam 
direction, cos*((9„), ensures the reconstructed neutrino 
is within the detector acceptance to reject events with 
significant radiation along the beam direction. The an- 
gle in the CM between the direction of the neutrino and 
the r, 0*„, is used to reject background events with a 
back-to-back track topology such as dimuon and Bhabha 
production. An event is accepted if it falls within a two- 
dimensional region defined with respect to m T and the 
e + //i + CM momentum, Events in this region are 
then used in a maximum likelihood fit to extract the sig- 
nal yield. 

The values of the selection criteria are shown in Ta- 
ble [J We optimize the selection sensitivity by defining a 
nominal signal box with a width of three standard devia- 
tions in the reconstructed m T andp*. The resolutions on 
m T and p*[ are approximately 10MeV/c 2 and 45MeV/c, 
respectively. The values of the selection criteria are cho- 
sen to maximize the discriminant S/ \[B where S is the 
number of MC signal events in the nominal signal box 
and B is the number of data events accepted outside this 



TABLE I: Selection criteria for the decay modes. The same 
criteria are used for the e + and p + lepton flavors except for 
E* ■ 

^miss • 



e + e — > It 
t~ 7r~ 7r + n~ v T 



— > fi t (e t 



?»iss (GeV) 0.015 
- 2 (GeV 2 /c 4 ) 



"mis: 

^ (GeV) 

COS*f0„) 



-0.9 



> 



m T (GeV/c 2 
Pi (GeV/c) 



1.6 
4.90 



3.23 

0.56 

0.20 

0.9 

0.015 

2.0 

5.32 



0.65 



-0.9 

1.6 
5.02 



4.55 (4.0) 

0.65 

0.15 

0.7 

0.090 

2.0 

5.32 



region but within 1.5 < m T < 2.2 GeV/ c 2 and the p\ 
boundaries given in Table HI We repeated the procedure 
using background MC within the nominal signal box in- 
stead of data and this produced consistent results. The 
signal MC reconstruction efficiencies and their statistical 
error after the application of these selection criteria are 
shown in Table HT1 

The backgrounds are dominated by e + e~ — * t + t~ 
decays where one r decays to an e + //i + plus neutri- 
nos and the other to either 7r~7r + 7r~^ r or n~v T . Light 
quark continuum processes are predicted to contribute 
significantly to e + e~ — > /i + t~(t~ — > 7r _ 7T + 7r _ ^ r ) only 



and events from e + e 



> fx n~~ are only present in 
(t~ — * n~v T ). Charm and BB back- 
grounds are eliminated by the track multiplicity and AE 
requirement and all other backgrounds are negligible. 

An extended unbinned maximum likelihood (ML) fit 
to the variables m T and p* is used to extract the total 
number of signal and background events separately for 
each mode. The likelihood function L is: 



N 



L = 



Nl 



(2) 



where rij is the yield of events of hypothesis j (signal or 
background) and TV is the number of events in the sam- 
ple. The individual background components comprise 
e + e~ — > t + t~, e + e~ — > fi + \i~ (7) and light quark con- 
tinuum decay modes. Vj{%i) is the corresponding prob- 
ability density function (PDF), evaluated with the vari- 
ables Xi — {m T ,p*} for the ith event. For the signal, 
we use double Crystal Ball functions [3] for both m T 
and p* . Due to correlations between m T and p* for non- 
signal events, we use a two-dimensional non-parametric 
PDF obtained from MC for the backgrounds [20( . In the 
maximum likelihood fit to the data, the parameters of the 
PDFs are fixed to the values determined from MC and 
only the signal and the background component yields are 
allowed to float. The statistical errors on the yields by 
the ML fit are roughly a factor of two smaller than those 
achievable with a simple counting experiment. 
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We check the robustness of the fitting procedure 
against variations in the signal size and background 
shape. We first fit the data outside the signal region with 
the MC background PDFs only, to determine their am- 
plitudes. Using these PDFs for the background, we gen- 
erate trial distributions including a Poisson-distributed 
number of simulated signal events, and perform the fit 
for each. We use 1000 trials at each of twenty values of 
the average signal yield between and 100 events, and 
find the fitted signal yield to be unbiased and the statis- 
tical uncertainty to be estimated correctly. Secondly, we 
generate a set of trial distributions in which the relative 
amplitudes of the simulated background components are 
changed, and confirm that this does not bias the fitted 
signal yield. 

As a validation check, we compare the predicted MC 
background levels and distributions of the variables from 
Table U to the data in the region outside the nominal sig- 
nal box and find that they are in agreement. We also 
extrapolate the fitted background PDFs from the region 
outside the nominal signal region into the nominal signal 
region and predict (measure) 193 ± 9 (202) and 143 ± 7 
(154) for e + e~ — > /i + t~(t~ — > 7r~7r + 7r~zy r ) and e + e~ — > 
/j + t~(t~ — > ir~is T ), respectively, and 112 ± 7 (128) and 



90 ± 6 (75) events for 



7T 7T + 7T V T ) 



and e + e~ 



(r — > 7r v T ), respectively, where the 



error is statistical only. The predicted and measured val- 
ues are consistent within the statistical errors. 

From the reconstructed MC efficiency, we can estimate 
the predicted number of background events and compare 



to the results of the ML fit. For e + e — ► r + 
predicted (ML fitted) background in the fitted region 
is 750 ± 43 (775 ± 19) and 494 ± 40 (385 ± 35) events 
for e + e~ — » fi + T~(r~ — > 7r~ 7r + 7r~ v t ) and e + e~ — > 
/i + T~(-r~ — ► k~v t ), respectively, and 414±41 (518±41) 
and 319 ± 45 (331 ± 18) events for e+e" -> e + T-(r- -> 
7r~ tt + 7r~ v T ) and e + e~ — > e + T~(r~ — > 7r~^ T ), respec- 
tively. The dimuon background to e + e~ — > /j + t~(t~ — > 
ti-v t ) is predicted (ML fitted) to be 114±38 (189±30). 
For the light continuum background, the MC predicts 
(ML fitted) 119±24 (129±40) and 19±9 (18±35) events 



/i + T (t 



7r 7T + 7r v T ) and 



e + T~(r _ — ► 7r~ tt + v T ) , respectively. The predicted 
and fitted values agree within errors. 

The main sources of systematic error on the signal 
yield come from uncertainties in the reconstruction, the 
r~ — > 7r~7r + 7r~z/ r decay mechanism and the fit pro- 
cedure. A relative systematic uncertainty of 0.8% per 
track, added linearly for all charged tracks in the event, 
is applied to account for differences in MC and data 
charged particle reconstruction. A relative systematic 
uncertainty of 1.0% per charged pion track and 1.3% per 
e + //i + track, added linearly for each charged track, is ap- 
plied to account for differences in MC and data particle 
identification efficiencies. 

A possible non-axial-vector decay mechanism for the 



decay r~ — > 7r~7r + 7r~^ r is not completely ruled out by 
current measurements 0]. To estimate this effect, the 
signal MC events were generated with 90% axial-vector 
and 10% phase-space decays and the difference in the re- 
construction efficiency compared to 100% a^(1260) de- 
cays applied as a systematic. This introduces a relative 
systematic uncertainty of 3.2%. 

The largest systematic error come from the variation 
of the PDF fit parameters within their fitted errors. The 
two-dimensional non-parametric background PDFs show 
small structures that depend on MC statistics and the 
value of the smoothing parameter used [io| . By vary- 
ing the smoothing parameter, using different functional 
forms and varying the fitted parameters within their un- 
certainties, we derive a systematic error of ~ 0.5 events. 
To investigate possible mismodeling of the detector ac- 
ceptance and response, we repeat the analysis with each 
selection criterion varied by the resolution on the cor- 
responding variable. All changes to the signal yield are 
smaller than the statistical error and we conservatively 
take the largest change in each case as a systematic un- 
certainty, which ranges from 2.5 to 4.4 events. The total 
systematic error is between 2.6 and 4.4 events and our 
final limit on the cross sections is dominated by the sta- 
tistical error which is of the order of 10 events. 

The m T and distributions for the modes are shown 
in Figure Q] and the projections are shown in Figures [2] 
and [3] The projection of the signal PDF is shown as the 
dashed line, the background PDFs as the dotted line and 
the total PDF as the solid line. The central value of the 



the cross section for e + e~ — > 1 + t~ is given by a — N/r]e£ 
where N is the number of signal events, r\ the signal re- 
construction efficiency and e is the r~ — > tt~tt + tt~is t or 
t~ — > n~v T branching fraction. The measurements are 
not statistically different from the null hypothesis and 
we obtain 90% CL upper limits by finding the maxi- 
mum number of signal events N such that the integral 
of the total likelihood function is 90% of the total in- 
tegral. From MC studies 12j, the total cross section 
of the process e+e~ — > fJ. + n~ at y/s = 10.58 GeV is 
a w = (1-13 ± 0.02) nb and we use this to calculate 90% 
CL upper limits on the ratio of the cross sections with 
respect to the dimuon cross section. The central values 
of the signal yields from the maximum likelihood fit and 
the upper limits on the cross sections and cross section 
ratios are given in Table ITT1 

We combine the r _ — > 7r~7r + 7r~^ r and r~ — > 7r~z/ r de- 
cays and calculate 90% CL upper limits on the cross sec- 
tions of a^r < 3.8 fb for e + e~ — > and a eT < 9.2 fb 
for e + e~ — > e + r~. The 90% CL upper limits on the ra- 
tio of the cross sections with respect to the dimuon cross 
section are calculated to be a^r/a^ < 3.4 x 10~ 6 and 



Cer/cW < 8.2 x 10 6 . For comparison with previous 
LEP results measured at > 92 GeV, the 95% CL up- 
per limits on the cross sections and ratio of cross sections 
are 4.6 fb and 4.0 x 10~ 6 for e+e" -> h+t~ and 10.1 fb 



FIG. 1: m T versus p* for reconstructed candidates for: a) 
e + e _ — > /i + T~(r~ — > n~ tt + tt~ v T )\ b) e + e~ — > ^i + t~(t _ — > 
tt~v t ); c) e + e~ — > e + T~(r~ — > ir~ it + -k~ v T ); and d) e + e~ — > 
e + T~(r~ — > 7r~^ T ). 




FIG. 2: Reconstructed distributions for e e~ — > /U r _ can- 
didates: a) m T and b) p£ for t~ — > 7r~7r + 7r~z/ T ; and c) m T 
and d) p* for r~ — > n~v T . The projection of the ML fit (solid 
line) hides the background component (dotted line) . The pro- 
jection of the few signal events is shown on the horizontal axis 
as a dashed line. The peaking dotted line shows the expected 
MC signal distribution at the 90% CL upper limit. 



and 8.9 x 10 6 for e + e — > e + r , respectively. 

In conclusion, we have performed the first search 
at a CM energy of y/s = 10.58 GeV of the lepton- 
flavor-violating production processes e + e~ — ► /i + r~ and 
e + e~ — > e + r~. No statistically significant signal events 
were observed in any of the decay modes. Upper limits 
have been placed on the cross sections and ratios of cross 
sections to the dimuon cross section to form limits on 
e + e~ — > /i + t~ and e + e _ — > e + r~. 

We are grateful for the excellent luminosity and ma- 
chine conditions provided by our PEP-II colleagues, and 



FIG. 3: Reconstructed distributions for e + e~ — » e + r~ can- 
didates: a) m T and b) p* for r~ — > 7r~7r + 7r~^ T ; and c) m T 
and d) p* for r~ — > ix~v T . The solid line is the projection of 
the ML fit, the dotted line is the background component and 
the dashed line is the signal component. The peaking dotted 
line shows the expected MC signal distribution at the 90% 
CL upper limit. 



TABLE II: Summary of the signal yields, cross sections and 
ratios of cross sections to dimuon cross section. The first 
uncertainty is the statistical error and the second systematic. 



e + e~ — ► /i + T~ 


T — > 7T 7T + 7T h> T 


T — > TT V T 


Total Events 


905 


575 


Signal Events 


-1.37 ±9.9 ±2.6 


1.9 ±10.1 ±4.4 


Signal Events (90% CL) < 19.2 


< 19.9 


MC Efficiency (%) 


18.5 ±0.2 


9.62 ±0.14 


cr^r (fb) 


-0.35 ±2.6 ±0.7 


0.85 ± 4.5 ± 2.0 


(7 fj,T (90% CL) 


< 4.9 fb 


< 8.9 fb 


a MT (95% CL) 


< 5.91 fb 


< 11.4 fb 


<W<7 MM (90% CL) 


< 4.3 x 10~ 6 


< 7.9 x 10" 6 


<W<r MM (95% CL) 


< 5.2 x 10~ 6 


< 10.1 x 10~ 6 




e + e~ — + e + r~ 


T — > 7T 7T + 7T V T 


r — > ty u T 


Total Events 


537 


332 


Signal Events 


15.9 ± 10.3 ±2.7 


10.7 ±8.8 ±2.7 


Signal Events (90% CL) < 32.3 


< 25.8 


MC Efficiency (%) 


11.73 ±0.15 


11.9 ±0.15 


<r eT (fb) 


6.5 ±4.2 ±1.1 


3.9 ±3.2 ±1.0 


a eT (90% CL) 


< 13.2 fb 


< 9.4fb 


a eT (95% CL) 


< 14.8 fb 


< 11.1 fb 


(Jer/O^ (90% CL) 


< 11.7 x 10~ 6 


< 8.4 x 10~ 6 


<7er/a^ (95% CL) 


< 13.1 x 10~ 6 


< 9.8 x 10" 6 
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We report on a search for the lepton-flavor-violating processes e + e~ — > /i + r~ and e + e~ — > e + r~. The data 
sample corresponds to an integrated luminosity of 211 fb -1 recorded by the BaBar experiment at the SLAC PEP-II 
asymmetric-energy B Factory at a center-of-mass energy of y/s — 10.58 GeV. We find no evidence for a signal and set 
the 90% confidence level upper limits on the cross sections to be < 3.8 fb and a eT < 9.2 fb. The ratio of the cross 
sections with respect to the dimuon cross section are measured to be a^rja^ < 3.4 x 10~ 6 and Uer/a^ < 8.2 x 10~ 6 . 



